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Spatially and temporally inhomogeneous evolution of one-dimensional vicious walkers with wall 
restriction is studied. We show that its continuum version is equivalent with a noncolliding system 
of stochastic processes called Brownian meanders. Here the Brownian meander is a temporally 
inhomogeneous process introduced by Yor as a transform of the Bessel process that is a motion of 
radial coordinate of the three-dimensional Brownian motion represented in the spherical coordinates. 
It is proved that the spatial distribution of vicious walkers with a wall at the origin can be described 
by the eigenvalue-statistics of Gaussian ensembles of Bogoliubov-deGennes Hamiltonians of the 
mean-field theory of superconductivity, which have the particle-hole symmetry. We report that the 
time evolution of the present stochastic process is fully characterized by the change of symmetry 
classes from the type C to the type CI in the nonstandard classes of random matrix theory of 
Altland and Zirnbauer. The relation between the noncolliding systems of the generalized meanders 
of Yor, which are associated with the even-dimensional Bessel processes, and the chiral random 
matrix theory is also clarified. 

PACS numbers: 05.40.-a, 02.50.Ey, 05.50.+q 



I. INTRODUCTION 

Unbalance between short-ranged properties of inter- 
actions among elements and long-ranged cooperative ef- 
fects realized in macroscopic levels is a significant fea- 
ture of systems far from equilibrium. Even in one di- 
mension the contact process, a model of infection of a 
contagious disease, exhibits a continuous phase transi- 
tion at a critical value A c of the infection rate A and in 
A > A c infected and healthy individuals establish coex- 
istence without detailed balance Q, @ . Boundary condi- 
tions locally imposed at the two edges in one-dimensional 
lattice play essential role to determine the bulk proper- 
ties in the asymmetric simple exclusion process, which 
can be regarded as a model of traffic flows in highways 
[J 0] • The purpose of the present paper is to propose 
one theoretical treatment of such emergence of long-range 
effects in simple (i.e. short-ranged) stochastic models us- 
ing vicious walker models originally introduced by Fisher 
for wetting and melting transitions [(|. The key-point is 
the symmetry of higher-dimensional space, in which the 
nonequilibrium many-body system is embedded. 

Consider N identical and independent simple and sym- 
metric random walks with initial positions x\ < x-i < 
■ ■ ■ < xn , where assumed to be even integers. One 
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of the fundamental quantities in the vicious walk prob- 
lem is the probability A/jv (t, x), x = (xi, • • • , xn), that 
all walkers remain the ordering of their positions up to 
time t; x\(s) < ^(s) < ■•• < %n(s) for all < s < t. 
In other words, it is the probability that they never col- 
lide with each other for a time period t. If two of them 
collide, then both are annihilated, since all walkers are 
vicious persons. It should be noted that this noncolliding 
condition seems to be very local and incidental, since any 
walker can enjoy free walking, while the relative distances 
from the nearest-neighbor walkers are greater than two 
units of lattice spacing. Fisher and Huse and Fisher 
Q derived the asymptotic form A/jv(i, x ) ~ t~^ N in large 
t for finite x and determined the exponent as 

^N = ^i 2 j=^N(N-l). (1) 

Interesting and important fact is that ipN is nonlinear in 
N expressing the long-ranged effect among vicious walk- 
ers, which is a result of accumulation of contact repul- 
sive interactions between nearest-neighbor walkers dur- 
ing the time interval [0, t]. Moreover, Q implies that 
the system possesses symmetry with respect to permu- 
tations of the walker positions. This hidden symme- 
try was clarified as explained follows. Huse and Fisher 
6, 8] mapped the enumeration problem of walks of the 
N particles from a set of positions xi < ■ ■ ■ < xn to 
Hi < • • • < TJn in time period t onto the diffusion prob- 
lem of a single particle in the N dimensional space with a 
set of wall restrictions (the phase space) from a position 
x = (xt, • • • , xn) to y — (yi, • • • , un) in time t. Assume 
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that p(t,y\x) denotes the transition probability density 
of a one-dimensional Brownian motion from x to y in 
time t, that is, p(t, y\x) = e~( y ~ x ' ' 2t /s/2nt, then by ex- 
ploiting the method of images, they derived the A^-body 
Green function of vicious walkers in the determinantal 
form [9| 

/(*,y|x)= i del ■, Jp(t,y k \ Xj )]. (2) 

1<j, k<N 

As discussed in 0, Q and explicitly shown in 0, (|2J) 
is found to be factorized into a product of symmetric 
part (the Schur function multiplied by Gaussian kernels) 
and the antisymmetric part (the product of differences of 
variables), and Af^(t,x) is obtained as an integration of 
PJl over y with restriction y\ <yi < • ■ ■ < yjv- 

Based on the above knowledge on the function 
A/jv(£, x), let us next consider the evolution of vicious 
walkers in time t, conditioned that they remain their 
ordering (i.e. noncolliding condition) up to a given fi- 
nite time T. Katori and Tanemura showed that this 
stochastic process was inhomogeneous both in space and 
time and a transition in the particle distribution was ob- 
served as time t goes on from to T. This transition 
is characterized by a symmetry change, which can be 
described not in the real one-dimensional space nor the 
Af-dimensional phase space but in the space of N x N 
hermitian matrices. That is, the problem was exactly 
mapped to the statistics of N real eigenvalues oi N x N 
hermitian random matrices in a time-dependent Gaus- 
sian ensemble. Due to the hermitian condition on N x N 
matrices with complex elements Hjk = H^ k + iHj k ,i 
V ' 1.1 < j,k < N, N(N + l)/2 variables in a set 
K = {Hf k : 1 < j < k < N} and N(N - l)/2 vari- 
ables in a set 1 = {H) k : 1 < j < k < N} are chosen as 
independent variables. These A^ 2 variables in total are 
assumed to be independently distributed following the 
Gaussian distributions with zero means. The variances 
of the variables in 1Z and X are proportional to ctr and 
<7i, respectively, both of which are functions of t. As the 
time t is approaching the final time T, the variance o\ 
decreases to zero and a transition from the ensemble of 
complex hermitian matrices (the Gaussian unitary en- 
semble, GUE) to that of real symmetric matrices (the 
Gaussian orthogonal ensemble, GOE) occurs. By inte- 
grating over the N' = N 2 — N other variables than the 
N eigenvalues, a transition of the eigenvalue-statistics 
from GUE class to GOE class is formulated 11] and it 
is indeed realized as the time evolution of positions of 
vicious walkers 0, 0, ^| . 

The above results suggest the possibility that vicious- 
walker-type problems with A^ walkers are generally 
mapped to some solvable problems in the spaces with 
appropriately higher dimensions N + N' , in which only 
the symmetries of the spaces should be considered and 
the interactions (restrictions) among the original walkers 
are resulted from integrating over the auxiliary A^' vari- 
ables. The symmetries of the higher dimensional spaces 
govern the macroscopic behaviors of the systems. The 
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FIG. 1: Vicious walkers with a wall. 



interacting particle systems far from equilibrium will be 
exactly solved, if we are able to find relevant symmetries, 
which are generally hidden in the original descriptions of 
the systems. 

Now we propose two kinds of problems of vicious walks, 
which will be solved in the present paper in order to 
demonstrate the above mentioned scheme for nonequilib- 
rium systems. We assume that all walkers are located in 
the positive region of position as < x% < x^ < ■ ■ ■ < xjy 
and put an absorbing wall at the origin (see Fig.l). The 
problems are (i) to determine the probability A/jv(i,x) 
that all walkers remain the ordering of their positions 
(noncolliding condition) with keeping apart from the wall 
up to time t, and (ii) to find out the time-dependent 
matrix model, whose eigenvalue-statistics realizes this 
stochastic process of vicious walkers with a wall, and clar- 
ify the hidden symmetry transition in the time evolution. 

The former problem (i) was already solved by Kratten- 
thaler et. al [lj| and the exponent governing the asymp- 
totic form in large t, A/jv(t, x) ~ t~'^ N , was determined 
as 

N 2 

4n = — ■ (3) 

So in this paper, we will start from their result and take 
the continuum limit of the model to solve the latter prob- 
lem (ii). In Section II, we construct a system of noncol- 
liding Brownian motions with a wall as a diffusion scaling 
limit of the corresponding vicious walker model. An im- 
portant result is that the A^-body Green function of the 
obtained system is also in the determinantal form @ for 
< x\ < ■■■ < xn,Q < yi < ■ ■ ■ < j/jv, if we replace 
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(a) (b) 

FIG. 2: (a) Three-dimensional Bessel process, (b) Method of 
images. 

P(t,y\x) by 

f( t , y\x) = -±= \ e - { y-^l 2t - e -fe+-) 2 /2*| . (4) 

This function becomes zero as y — > 0, representing the 
effect of the absorbing wall at 0. There are two dis- 
tinct ways to derive ijlfl. (a) Consider a Brownian mo- 
tion started from the origin in a space with dimensions 
d > 2. We adopt the d-dimensional spherical coordinate 
r = (r, 6±, ■ ■ ■ , 6d-i) to represent the motion. In partic- 
ular, we can trace the radial coordinate (the modulus of 
the Brownian motion) r = r(t) as shown in Fig. 2 (a) for 
d = 3. Since the transition probability density of r is gen- 
erally described using the Bessel function, such a stochas- 
tic process of r is called the Bessel process 0, 0, ■ 
If we multiply the transition probability density of the 
three-dimensional Bessel process by x/y, then Q is ob- 
tained, (b) For a real path of the Brownian motion from 
(x, 0) to (y, t) in a spatio-temporal plane, we consider an 
imaginary path from (—a;, 0) to (y, 0), where x, y > (see 
Fig. 2 (b)). As an analogy of electrostatic problem, we 
subtract the transition probability density of the imagi- 
nary paths from that of the real paths to obtain 0} (the 
method of images). 

Yor studied a temporally inhomogeneous process called 
the Brownian meander, which is obtained as a trans- 
form of the three-dimensional Bessel process used in the 
above derivation (a) of (0J. He also introduced the d- 
dimensional generalized meanders as the transform of the 
d-dimensional Bessel processes E3- I n Section III, we 
will give a general theory of noncolliding N walkers con- 
structed as a conditioned system of d-dimensional gener- 
alized meanders. As a special case, it provides a proof 
of that the noncolliding system of Brownian particles in 
the presence of a wall is equivalent with the noncolliding 
system of the Brownian meanders. This is the complete 
generalization to many particle systems with arbitrary 
number of particles N of the fact that the single-particle 
Green function Q of a Brownian motion with wall re- 
striction at the origin is proportional to the transition 



probability density of a single three-dimensional Bessel 
process. A key point of our proof for this new result is 
the proper transform from the iV-body Green function 
@ to the transition probability density by multiplying 
an appropriate ratio of A/jv's (see Eqs. (JSJ and lf^U|l be- 
low). Moreover, our general argument gives that, if we 
consider the problem (ii) for the c?-dimensional Bessel 
processes and generalized meanders in the case of even d, 
it is solved using the Gaussian matrix theory with chiral 
symmetries, which is relevant for the physics of QCD at 
low energies El El USES- Since Brownian meander is 
made from the three-dimensional Bessel process as men- 
tioned above, it is concluded that the present problem 
with Q is not related with chiral symmetry of matrices 
and that matrix models in the different symmetry classes 
should be considered. 

The latter derivation (b) of Q gave us a hint to find 
out the true symmetries, which govern the distribution of 
vicious walkers with a wall; particle-hole symmetry. The 
particle-hole symmetry is important in the BCS theory of 
superconductivity. In particular, its microscopic mean- 
field treatment ignores any local interactions among par- 
ticles and holes but consider this symmetry with the 
so-called Bogoliubov-deGennes (BdG for short) Hamilto- 
nian. Random matrix theory of the BdG-type Hamilto- 
nians was introduced and developed by Altland and Zirn- 
bauer in order to describe the energy-level statistics and 
transport properties in a metallic quantum dot in contact 
with a superconductor in a magnetic field |22l |23| . They 
studied the Gaussian ensembles of the BdG Hamiltonian 
in the form 

W = C/t (o -u) U ' W = dia s( w i' ' ' ' '^n), 

where U is an appropriate unitary matrix. As we will give 
a brief summary of their results in Section IV, Altland 
and Zirnbauer discovered four new classes of eigenvalue- 
statistics in addition to the previously known three classic 
Wigner-Dyson classes (two of them are GUE and GOE, 
whose relation with the vicious walks without wall was 
reported in [10|) and three chiral symmetry classes (two 
of them will be argued in Section III associated with the 
noncolliding systems of generalized meanders constructed 
from the even-dimensional Bessel processes). In the four 
nonstandard symmetry classes, denoted by C, CI, D and 
Dili in Cartan's notation [2^, the classes C and CI are 
relevant in the present vicious walk problem. They have 
the probability density functions of nonnegative eigenval- 
ues in the form 

N 

p^W)cxe-M 2 /^ Yl | Ufc - u /IlN a . 

l<j<k<N 1=1 

(5) 

where \lo\ 2 — YljLi an d the indices a and are spec- 
ified as 

a = 2, [3 = 2, for class C, 
a = 1, fi = 1, for class CI 
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respectively. We will show in Section IV that the tran- 
sition of distribution of vicious walker positions with a 
wall is described by the symmetry change from the class 
C to the class CI of the BdG Hamiltonians. This fact was 
already reported by Nagao [24|, but in this paper com- 
plementary new results will be given. In an earlier paper 
[l3| , the transition from GUE to GOE realized in the vi- 
cious walks without wall was characterized by the graph- 
ical expansions with time-dependent coefficients for the 
moments of walkers. In Section V, we will introduce the 
Mobius graph expansions for the moments of the vicious 
walkers with a wall. Moreover, using an exact results 
of dynamical correlations by Nagao [24| , closed formulae 
for the moments will be given. Such graphical expansions 
will be first reported in detail in the present paper for the 
nonstandard symmetry classes C and CI of Altland and 
Zirnbauer. Concluding remarks are given in Section VI. 



II. VICIOUS WALK WITH A WALL AND ITS 
DIFFUSION SCALING LIMIT 

A. Determinantal formula 

First we consider the N independent simple and 
symmetric random walks on an integer lattice Z = 
{• • • , —2, —1, 0, 1, 2, • • •} starting from the sites {xj}, j — 
1, 2, • • • , N, and denote the position of j-th walker at time 
n = 0, 1, 2, • • ■ by Xj (n). Assume that the initial positions 
are all distinct nonnegative even integers and ordered as 
< xi < X2 < ■ ■ ■ < xn- Then we impose the noncollid- 
ing condition up to a given time m > 0, 



Xi(n) < X2(n) < ■■■ < xjv(n), 



n 



1,2,- 



, TO. 



Such conditional walks are called vicious walks up to time 
to 0. Here we impose further restriction on the walks 
as 

Xj(n) > 0, j = 1, 2, • • • ,N, re = 1, 2, • ■ • ,m. 

In other words, there is a wall at the origin and all the 
walkers are conditioned never to collide with each other 
nor to collide with the wall during the time interval < 
n < to. 

Let N N (m,y\x), x = • • • , x N ),y = {y x , ■ ■ ■ ,y N ), 
be the total number of the vicious walks with wall re- 
striction, in which the N walkers start from the po- 
sitions Xj , j — 1,2, • • • , iV, and arrive at the positions 
yj , j = 1, 2, • • • , N, at time to. Krattenthaler et al. gave 
the determinantal formula to this number as 141 



N N (m,y\x) 



det 

l<j,k<N 



m 

m+Xj-y k 



rn 

m+Xj+yk + 2 
2 



Suppose that all random walks start from given initial 
positions x. Since the total number of walks is 2 mN , the 
probability that they are vicious walks with wall restric- 
tion and end up with positions y is N]y(m, y|x)/2 mAr . 



B. Diffusion scaling limit 

In order to take the continuum limit of the vicious 
walks to derive the system of noncolliding Brownian mo- 
tions, we introduce a function 4>l(x) — 2[Lx/2], for 
L > 0,x £ R (the set of all real numbers), where [z] 
denotes the largest integer not greater than z, and let 
0i(x) = {4>l{x\), ■ ■ ■ , <Pl(xn))- By Stirling formula, we 
can take the diffusion scaling limit as 

(L\ N 1 ~ 

l^oo {2) ^N NN ^(t),My)\<t>L^)) = Mt,y|x) 

for x, y with < Xi < ■ ■ ■ < xn, < yi < ■ ■ ■ < un, 
where 

f N (t,y\x) = det fflt,y k \xj)] (6) 

l<J,k<N 

with J3J . Let T > and consider a system of N Brownian 
motions conditioned never to collide with each other nor 
to collide with the wall at x = in [0, T]. Set 



M N {t,x) 



dyfN(t,y\x) 



(7) 



0<Vi<-<Vn 



where dy = ri, =1 dyj. Then the transition probability 
density from the state < Xx < ■ • • < xn at time s to 
the state < y\ < ■ ■ ■ < j/jv at time f(> s) of such a 
system is given by 

- / 1 \ fN(t~ s,y\x)Jv N (T- t,y) 
g N , T (s, x; t, y) = T — r , (8) 

since the numerator is the probability that we have the 
noncolliding (with each other and with a wall) Brownian 
paths from x at time s to y at time t and these paths keep 
noncolliding from the time t up to time T as well, and the 
denominator is the probability that the Brownian paths 
are noncolliding all during time interval [s,T]. 
It is useful to rewrite © as 

f N (t, y|x) - (2 7 ri)- iV / 2 S p ?(y) (W', • • • , e*"/*) 

N 

2tJJ( e ^/i_ e -^/i) 
3=1 

JJ /( e »i/t_ e *i/t)( e (»«+**)/t _i)J 

i<j<N 

N -N+l 

{n^ 7 *} . ( 9 ) 



x e 



l<i<j<N 
N 



3 = 1 



where £(y) = (^i(y), ■ • ■ ,^(y)) with ^-(y) = y N -j+i 
(N-j + l),j = l,2,--.,N, 



sp ? (zi, ■■■ ,z N ) 



det(z. 



Cj + N-j+l -{Hj+N-j+l) 



det(z. 



N-j+l _ -(N-j + l) 
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and |x| 2 = J2j=i x2 j- Remark that sp A (zi, • • • , zn) is the 
character of the irreducible representation corresponding 
to a partition A of the symplectic Lie algebra (see, for 
example, Lectures 6 and 24 in |25j). Since we know the 
formula 

£2 _ e2 n 

sp«(v--,i)= n ^An— 

l<i<j<N 3 1 j=l 3 

with ij = f 3 - + N - j + 1, m,j = N - j + 1 and the 
integral ((17.6.6) on page 354 in ptjjl 

. AT 
l<i<j<Af j=l 

r(i + j 7 )r(a + 7 (j-i)) 
JUL r(i + 7 ) 

we have the asymptotic 

Af N {t,x) = ^-h N (x/Vi) (l + 0(x/Vt)) (10) 

C N \ ) 

for |x|M - 0, where c^v = (tt/2)^ 2 nJLi r(2j)/r(j) 
with the Gamma function L(z) and 

AT 

Mx)= n (4-^ 2 )n^ ( n ) 

l<j<k<N 1=1 

Substituting eqs.® and lflH|l into eq.JHJ), we find 

g N , T (0,0;t,y)=c N T N2 /H- N ( 2N +W 2 

xe-^ 2 / 2t h N (y)rt N (T~t,y), (12) 

where cjv = 1/ Oj=i ^(j)* ^ ne transition probability 
densities (jHJl and (|12|l define the N noncolliding Brownian 
motions with wall restriction in time interval (0, T] |27j. 

It should be noted that A/iv(i,x) is the noncolliding 
probability of Brownian motions with wall restriction and 
(|10fl gives the power-law behavior A/Ar(i, x) ~ j-V'n i n 
large t for finite x with the critical exponent Q. 



where d N = (2/tt) jv / 2 / n£Li T(2j). That is, we have the 
identity 

p N (0,O;t,y) = Nlp™ G (y;t) 

for < y\ < ■ ■ ■ < un, where 2 G ( W ) <y2 ) 1S * ne proba- 
bility density function (J5J) of nonnegative eigenvalues of 
the BdG Hamiltonian in the class C (a = (3 = 2). On 
the other hand, if we set t = T in JHJl, for A/jv(0, y) = 1, 
we have 

^,T(0,O;T,y) =?^- w ( JV + 1 )/ 2 e-l"l 2 / 2i ^( y ), (15) 

which implies the identity 

g NtT (0,0;T,y) = Nlpff(y;t). 

That is, at t = T, 5^^(0,0; T,y) is identified with the 
probability density function J5J of nonnegative eigenval- 
ues of the BdG Hamiltonian in the class CI (a = (3 = 1). 

The above results mean the following facts. If we con- 
sider the N noncolliding Brownian motions with wall re- 
striction up to a finite time T > 0, in which all particles 
start from the origin, as the ratio t/T —* 0, the distribu- 
tion of particle positions is asymptotically described by 
the eigenvalue-statistics of the BdG Hamiltonian in the 
class C. On the other hand, at the final time t = T, 
it can be identified with the eigenvalue-statistics of the 
BdG Hamiltonian in the class CI. There occurs, thus, a 
transition from the class C distribution to the class CI 
distribution as time t goes on from to T in our stochas- 
tic process. 

The essential difference between (|14f) and <|15[) is found 
in the exponents of the factors /iAf(y) /3 such that (3 = 2 
in H4|) and (3 = 1 in Q15J1. This factor expresses strong 
repulsive interactions among particles and between the 
wall and each particle, in which the larger exponent (3 
gives stronger repulsion for short distances. At the very 
early stage of the process, t/T <C 1, the repulsion may be 
strong, since the noncolliding condition will be imposed 
for a long time period up to time T in the future. As the 
time t goes on, the repulsion strength is decreasing as is 
the remaining time until T, and attains its minimum at 
t = T. 



C. Transition from class C to class CI 

By IjlOp. the T — > oo limit of JSJ is determined and 
simply given as 

p N (0,x;t,y) = lim g NtT (0,x;t,y) 

T — ►oo 

= |^M*,y| x)) (13) 

where we have set s = and used I jllj l. Moreover, we 
can take the x — > limit of l|13(l to obtain 

p N (0, 0; t, y) = c< N t- N ^ N +V/ 2 e-M 2 / 2t h N (y)\ (14) 



D. Stochastic differential equations 

As explained in [T(il | in the case without wall restriction, 
the positions x(t) = (xi(t), ■ ■ ■ ,X]y{t)) of the N noncol- 
liding Brownian motions solve the stochastic differential 
equations in a modified type of Dyson's Brownian motion 
model |2(| [2^]. In the present case with wall restriction, 
we have 

dxj(t) = E?(x(t))dt + dBj(t), 
for < xi < ■ ■ ■ < x N , < t < T, where 

Ej(x) = -^-lnff N (T-t;xL), 



6 



and {Bj(t)}j =1 are N independent standard Brownian 
motions 

B,(0)=0, (B j (t))=Q, 

((B 3 (t) - B 3 (s)) (B k (t) - B k (s))) = \t- s\S jk , 

for any t, s > 0, j, k = 1, 2, • • • , N. In particular, in the 
limit T — > oo, 1(10(1 with (|llf> gives the equations 



dxj(t) = dBj(t) + 



E 

l<k<N,k^j 



Xj(t) 



-dt 



xj (t) - x k (t) Xj(t) +x k (t) 



dt (16) 



for 1 < j < N j27| ■ In the stochastic differential equation 
for the position of j-th particle, the drift terms dt/xj(t) 
and dt/(xj(t) — x k {t)) represent repulsive force from the 
wall at the origin and that from the fc-th particle, respec- 
tively. In addition to them, there are the terms in the 
form dt/(xj(t) + x k (t)),l < k < N,k 7^ j, which can 
be interpreted as effective repulsive forces from the mir- 
ror images of other particles located at —x k (t), 1 < k < 
N,k^j. 



III. NONCOLLIDING MEANDERS AND 
CHIRAL RANDOM MATRIX THEORY 

A. Definitions of elementary processes 

Consider a diffusion equation in dimension d > 2 



1 d d 2 
-«(t,y|x) = -£_, u (t,y|x) 



dt ' 2^dy 2 

with the initial condition u(0, y|x) = <5(x — y). Wc 
use the spherical coordinates x = [x, 6\, ■ ■ ■ , 6><j-i), y = 
(?/, ipi, • • • , tpd—i) an d integrate over all the angular vari- 
ables to obtain a differential equation for the radial co- 
ordinate (the modulus) 



9 -, , N 1 
di U{tMx) = 2 



dy 2 



d- 1 d 
y dy 



u(t,y\x) 



The unique solution of this equation satisfying the initial 
condition u(0, y\x)y d ~ 1 dy = 5(x — y)dy for x > is given 
as 

fift *> = ^r , ' w,,2 Mf)> 

where v = (d — 2)/2 and I v {z) is the modified Bessel 
function 



m = E 



(z/2) 



2n+f 



n\T(v + n+ 1)' 



n=0 



If we set p^(t, y\x) = u(t,y\x)y d 1 , then it is normal- 
ized as J °° p( v \t,y\x)dy — 1 for any x > 0. We define 



p^\t,y\Q) by the x -> limit of p("\t,y\x). Then we 
have 



P ^(t,y\x) = y^h-(^v 2 V*%(^),x>0,y 



P (v) (t,y\0) = 7T„ 



..2v+l 



2 v V(v+l)t v + l 



, y> 0. 



> 0, 



(17) 



The d = 2(y + 1) dimensional Bessel process is defined 
so that its transition probability density is given by l|17|) 

EEG3. 

For < u < T, w > 0, n € [0, 2(v + 1)), we consider 



i^'^ (u, w) 



dz p {v) (T -u,z\w)z- K . 



(18) 



That is, we multiply a weight z~ K at the final time T, so 
that, as the arrival position z is nearer to the origin, the 
path is more enhanced. Then the transition probability 
density from x at time s to y at time t with such bias at 
time T, < s < t < T, is given by 



pp K '(s,x;t,y) 



1 



(s,x) 



p^(t-s,y\x)hp K) (t,y) 
(19) 

for x, y > 0. This bias makes the process defined by 119(1 
be temporally inhomogeneous, and Yor called it the gen- 
eralized meander indexed (y, k). In particular, when v = 
1/2 and K = 1, the process is called the Brownian mean- 
der 01 • The transformation from (|17fl to Pj! (s, x; t, y) 
by ((19(1 is a generalization of the ^.-transform of Doob 
[29|. Remark that, if we replace T by u in 

O, then 

hu' K \u, w) — w~ K , since p("\0, z\w) = S(z — w). In the 
case v = 1/2, k = 1 and s = 0, ((19(1 becomes, by this 
replacement, p^ 1 / 2 ^, j/|a;) x x/y, which is equal to J2J for 
A/2( z ) = -\/2/ (7rz) sinhz. This is the derivation (a) of 
Pi mentioned in Section I. 



B. Non-colliding systems and rectangular random 
matrices 

Now we consider a system of N generalized meanders 
conditioned that they never collide with each other for 
a time interval (0,T],T > 0. Using the determinantal 
formula in J2J) and following the same way as © and 
((19(1 . the transition probability density is given by 

sk T '(*,x;t,y) = ' : (2°) 

■^Ar.r i T ~ s > x ) 

for < s < t < T,0 < 11 < ••• < x N ,0 < yi < ■■■ < y N , 
where 



f { N,T (s> x ;*>y) = ,det 

l<j,fe< AT 



pp K) {s,Xj]t,y k ) , (21) 



and 



A^ T K) (T-i,x) = 



0<yi<-<j/jv 



dyf^\T-t,x,T,y), 
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Since f^'tp (s, x; t, y) is temporally homogeneous and 

independent of T, we will write it as (t — s,y|x). 
Moreover, note that 



Pn,t (s,x;t,y) = , ^- -/^(t - s,y|x)/i^' re ''(t,y), 

/ij.' (s,x) 

(22) 

where /i^ (t,x) = IljLi an d that 

hp K) (T,x) = n^Li a^7 K - Then (j23> can be written as 



1 



9n t ( s > x ;*>y) 



/^- S ,y|x)A/^(T-i,y) (23) 



■(">«)/ 



^ K) (T- S ,x) 



with 



<" K) (t,x) 



AT 



0<ki<"-<8/jv 



dy/W(t,yWn 



3/; 



(24) 



i=i 



The important point is that we can confirm that 123(1 with 
v = 1/2 and k = 1 is equal to the transition probability 
density (JSJ of the noncolliding Brownian motions with a 
wall. In other words, we found the equivalence between 
the noncolliding Brownian motions with a wall and the 
noncolliding Brownian meanders. If we set TV = 1 in i|24|) 

with v = 1/2, k = 1, we have N^ 12 '^ (t, x) = 1/x. Then 
(|2^|l with s = is reduced to the equality p(t, y\x) — 
pt 1 / 2 ) (t, y|x) x cc/y, which is the statement given as the 
derivation (a) of (J3J in Section I. 

We then consider the limit |x| — > to define the non- 
colliding generalized meanders all started from the origin 
at the initial time s = for nonnegative integers v (i.e. 
even d) . Consider an arbitrary Ni x N 2 complex matrices 
with Ni > N 2 . We denote by M(Nj_ N 2 ; C) the space 
of all such matrices. It is known |30j | that any matrix 



A 6 M(Ni,N 2 ; C) can be expressed by 
A = U^AV, 



where U and V are unitary matrices with sizes Ni and 
N 2 , respectively, and A is the Ni x N 2 matrix in the form 



A 



with A = diag(ai, a 2 , ■ ■ ■ , ajv 2 ), (26) 



and where aj > 0,1 < j < N 2 . The matrices (U,V) 
parameterize the coset space U(N±) x U(N 2 )/[U(l)] N2 , 
where [U(1)] N2 is the diagonal subgroup of U(N 2 ), and 
thus the (U, A, V) can be regarded as "spherical coor- 
dinates" in the space M(Ni, N 2 ;C). It should be re- 
marked that {eti, ■ • • ,ajv 2 } are not eigenvalues of A; they 
will be referred to as "radial coordinates" . The following 
integral formula proved in [2(| is useful. Let dfJ,(U, V) 
be the Haar measure of U(N{) x U(N 2 )/[U{1)] N2 . For 
A, Be M(N!,N 2 ;C), set A = U A A A V A , B = UlA B V B , 
where U A ,U B G U(Ni), V A ,V B e L r (AT 2 )/[C/(l)]^ ) 



A, = 



Aa 




Ar = 



A B 




with A A = diag(ai,---,ajv 2 ), A B = diag(bi, • • • , b N2 ), 
aj > 0,bj > 0,1 < j < N 2 . Then for an arbitrary 
constant <r, 

j dn(U A ,V A ) exp(-^tr{(A-B^(A-B)}\ 
deti< jik <N 2 exp (-^^J In!-n 2 



N 2 

3=1 l<j<k<N 2 

Remark that this integral formula can be regarded as a 
version of the Harish- Chandra (Itzykson-Zuber) formula 

Using this integral formula, l(21|l with a nonnegative 
integer v and n = is written as 

/IrWix) 



n ^-4)Uyr +i n 



l<j<k<N 



3 = 1 



l<j<k<N 



x J dn(U Y , V Y ) exp (-^tr{(A - Y)\X - Y)f) . 

Since tr{(X - Y)\X - Y)} -> trY^Y = |y| 2 as |x| -> 0, 
we have 



lim 



^ II ($-4) 

l<j<k<N 



N 



« rw +i n w-«2)e" bri /M - ( 2? ) 

j=l l<j<k<N 

Above argument proves the following result. Let v be a 
nonnegative integer and < k < 2(v + 1). The limit 



(25) l x l of g ( N y (0, x, t, y) is given by 



9n',t (°>°;*>y) = c e" 

TV 

x Lhr +1 II (^-^^(r-t.y), (28) 



= - P -|yl 2 /2t 



i=i 



l<j<k<N 



where c is a normalization constant determined by 

Jo<vx<---<vN < ty 9 N,T ( ' ;*>y) = L The ^ noncolliding 
generalized meanders all started from the origin at time 
is defined by the transition probability density (|28|) . 



C. Chiral Gaussian ensembles and transition of 
chiral symmetries 

For the space M(N\,N 2 ;C) of all N\ x N 2 complex 
matrices, we introduce the integration measure dv(A) = 
IlfliUtidAf k dA] k for A = (A jk ) e M(N U N 2 ;C) 
with Ajk = Aj" k + iA l j k ,i = y— T. The chiral Gaussian 
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unitary ensemble (chGUE for short) with variance a 2 is 
the ensemble of matrices A £ A4(Nx, N 2 ; C) with the 
probability measure 

dfi chGUE {A;a 2 ) cx exp (--^trA^ dv{A). 

For A £ M(Ni,N 2 ;C) with the polar coordinates lf2"5|) 
and (|26[) . we can show that the probability density func- 
tion of the "radial coordinates" a = (ai, • • • , ajv 2 ) of 
A £ M.(Ni, N 2 ; C) in chGUE with variance a 2 is given 
as 



chGUE 



(a; a 2 ) 

JV 2 



-|a|V2 CT ' 



2(N 1 -N 2 ) + l 



n <« 



l<j<fc<JV 2 



Next we set M{N 1 ,N 2 ; R) be the space of all Ni x N 2 
real matrices for Nx > N 2 . The chiral Gaussian orthogo- 
nal ensemble (chGOE) with variance a 2 is the ensemble 
of matrices B £ A4(Ni, N 2 ; R) with the probability mea- 
sure 

dn chGOE {B;a 2 ) cx exp {^^^^B^j dv'{B) 



with dv'(B) = ]lfc=i dB ]k- The probability density 

function of the radial coordinates b = (61 , • • • , bjsr 2 ) is 
given in the form 



jN 2 



N 2 

p chGOE (b;a 2 ) oc e -l b l 2 /2- 2 Y[ b N.- 



N 2 



n 

l<j<fe<AT 2 



Then we consider the distribution of the sum of two 
rectangular matrices C = A + B, in which A and B 
are chosen from chGUE and chGOE, respectively. The 
distribution function of C is the convolution of those of 
chGUE and chGOE. Consider the ensemble of matrices 
C £ A4(Ni, N 2 ; C), in which the probability measure is 
given as 



^chGUE/GOE (c . CT 2 j(J 



d/i chGUE (C - B- a\)d^ GOE (B- a 2 ). 

BeM(N 1 ,N 2 ;R) 

We denote the probability density function of the radial 
coordinates c = (ci, • • • , cjv 2 ) of matrix C in this ensem- 
ble by p^ue/go^cjct 2 ,^; A^, AT 2 ). 

Comparing the above definitions and l|28l) . we can 
prove the following equality for nonnegative integers is, 



3 ^ +1) (0,O;t,y) 
N\ p chGUE / GOE (y;i(l 



T J 1 T 



N + v,N 



where < y\ < • • • < un- It implies that, if v = 
0, 1, 2, • • •, the time-evolution of the noncolliding gener- 
alized meanders indexed (V, v + 1) is represented by the 
transition of the eigenvalue-statistics from the chGUE 
class to the chGOE class. 



IV. GAUSSIAN ENSEMBLES OF 
BOGOLIUBOV-DEGENNES RANDOM 
MATRICES 

Since we have found that the noncolliding Brownian 
motion with a wall is equivalent with the noncolliding 
system of Brownian meanders with indices v = 1/2 
and k = 1, it does not belong to the chiral symmetry 
classes discussed in the previous section. We have to con- 
sider the ensembles of hermitian matrices in the form of 
the Bogoliubov-deGennes Hamiltonian for the mean-field 
theory of superconductivity. 

As shown in H-C, the time evolution of spatial distri- 
bution of particles in the present system can be regarded 
as a transition of the eigenvalue-statistics of the BdG 
random Hamiltonian in the class C to the class CI. The 
former statistics is characterized by the exponent (3 = 2 
of the repulsive factor h^iy) and the latter by j3 = 1. 
It does not imply, however, that the functional form of 
the distribution is maintained in the form JSJ) with a = [3 
and only the exponent [3 changes continuously as the time 
passes. It this section, using a version of Harish-Chandra 
(Itzykson-Zuber) integral formula over unitary group, we 
show the fact that the time evolution of the present pro- 
cess is described by a two-matrix model coupling random 
matrices, one of which is chosen from a Gaussian ensem- 
ble of the BdG Hamiltonian matrices of class C and other 
of which is from that of class CI. There the time depen- 
dence of variances of these two ensembles is different from 
each other. 



A. Hermitian and real symmetric matrices with 
particle-hole symmetry 

We consider the space of the hermitian matrices spec- 
ified by the following: 



.Vl"'"T2.Y:Cl = <iW=( * t _ b aT 



a is an N x N hermitian matrix, 

and b is an N x complex symmetric matrix X. 

Since the dimension of the space of a is N 2 and that of b is 
N(N+ 1), the dimension of X BdG (2iV; C) is N(2N+1). 
Define 



= 1 -I 



(29) 



where In is the N x N unit matrix. Then H £ 
M BdG (2N; C) has the following symmetry |H|2| 



H = -CWC 



T n -l 



(30) 



Assume that ipj is the 2N x 1 eigenvector of H 
with an eigenvalue 10 j\ Titpj — tOjipj. Then by i|3U[l . 
—C'H T C~ 1 (pj = Uj<Pj- Take the complex conjugate of 
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the both sides and use the hermiticity of Ti and the fact 
C- 1 = -C, we have H(Cip*) = -uj-{Cud*). That means 
C<p* is the eigenvector of H with the eigenvalue — ujj . As- 
sume that oj\, ll>2, ■ ■ ■ , ujn be the nonnegative eigenvalues 
of H, then other eigenvalues are given by — u>i, ■ ■ ■ , — ljn- 
Therefore, if we set 



u 



then 



HU^ = C/U, with A = 



Ui 



u 2 
-u* 



(31) 





— UJ 



We assume that U is unitary. Then we can see that ill 
satisfies the relation C = (iU) T C(iU). The set of such 
2N x 27V unitary matrices is called the symplectic group 
Sp(27V; C) [H |2l , whose dimension is £ = 2N 2 . 

The above consideration is summarized as follows; any 
H £ M BdG (2N; C) can be diagonalized as 



unu^ = a = 



uj 
-uj 



iU e Sp(2iV;C), 



where uj = diag(oji, • • • 
then consider the map 



,uj n ), ojj > 0,1 < j < N. We 



{u,U) = \ (wj)i<j< N ,p= (Pn)i<^<i^, (32) 



where p denotes the ^-dimensional vector, whose ele- 
ments are the independent variables of U. We have the 
Jacobian of this map as |23j . 



Jac(ip) = 



f an on an on 

\OUJi OUJpj Opi Opi 



N 



= c (p)iH n K 2 -^ 2 ) 2 > ( 33 ) 

j=l l<j<k<N 

where C(p) is a function independent of the eigenvalues 

UJ. 

Next we consider the set 



j\A (2N; R) 



ft 



a and & are N x N real symmetric matrices} . 

In this case, since the dimensions of the spaces of a and 
b are both N(N + l)/2, the dimension of M BdG (2N; R) 
is N(N + 1). We can see that any H £ M BdG (2N;R) 
can be diagonalized as 



U'HU' T = A = 



uj 
-uj 



iU' £ Sp(2AT;iR), 



where uj = diag(wi, • • • ,ujn), > 0, 1 < j < N. Here 
Sp(2A; iR) is the symplectic group of 2A~ x 2N matrices 
whose elements are pure imaginary. The map 



H 



°-+ (u,U') = ((wj)i<j<N,v' = (p M )i<p<i>), (34) 



with £' ^dimension of Sp(2N;R) = N 2 , is considered. 
Here p' denotes the ^'-dimensional vector with the ele- 
ments of the independent variables of U' . The Jacobian 
of this map is determined as [2j| 

JV 

Jac(^) = C"(p')ni^l II 1^-^*1. (35) 

j = l l<j<k<N 

where C'(p') is a function independent of the eigenvalues 
uj. The important point is that Jac(</?) and Jac(<//) are 
proportional to |/ijv(w)| 2 and |/ijv(w)|, respectively. 

B. Gaussian ensembles and eigenvalue distributions 

Altland and Zirnbauer introduced the Gaussian ensem- 
bles of the BdG Hamiltonians in M BdG (2N;C) and in 
M BdG (2N; R), in which the probability measures are m 
the form 

H N (H;a 2 )dH oc exp (-^2 tr ^ 2 ) (36) 

with variance 2cr 2 . For H £ M BdG (2N; C), we write 
the complex variables as ajk — a^ k + ia l j k ,bjk — b^ k + 
ib l - k with i = V^T, af k ,Oj k ,bf k , b l jk £ R, and choose the 
independent variables as {af k , bf k , b l - k : 1 < j < k < 
N} U {a] k : 1 < j < k < N}. Since 



trH 2 = \ H it 



N 



= 2E{(4) 2 + (^) 2 + (^) 2 } 

+ 4 £ {(af k f + {a) k f + (bf k f + (b) k f} 

l<j<k<N 

the probability measure l|36|l is rewritten as 
fi N (H;cr 2 ) 

e -(af 3 f/2^ e -(blf/2* 2 e -(b) 3 f/2** 



3=1 



V2^ 



V2 



7T(7 



V27> 



n 

l<j<k<N 

>< n 

l<j<k<N 



(37) 



'no 



with the integration measure 

N 

dH=l[daf j dbf j db I jj [] da%da) k db} k db) k . (38) 

3=1 l<j<k<N 

The probability measure H36|l is given for H £ 
j\4 BdG (2N; R), by setting all the imaginary parts 

X 3'fc' 6 3'J 



{ah k ,b l jk } be zeros in (J2U) and fSHft. 
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On the other hand, following the maps (|32() and (|34|l . 

the integration measures are transformed as 

dH oc 3a,c(ip)dui x dU oc h N {ui) 2 du)dU (39) 

for H E M BdG (2N; C), %U G Sp(2iV; C), and 

oc Jac(co')dw x d[/' oc \h N (uj)\dujdU' (40) 

for ft e X BdG (2iV;R), zf7' S Sp(2iV;iR), respec- 
tively. Since trft 2 — 2|w| 2 = 2^^ =1 wJ, integrating 
over the spaces Sp(2^V; C) and Sp(2iV; iR) gives the dis- 
tribution functions of the nonnegative eigenvalues lo = 



(u>x, ■ ■ ■ ,u!pf),u>i > 0, in the form of JSJ) with the indices 
a — (3 = 2 (class C) and with a = f3 = 1 (class CI), 
respectively. 



C. Harish-Chandra integral formula 

The transition probability density (|12l) from the state 
to the state y in time t is written as follows using J7J), 
where the proportional constants are independent of the 
stochastic variables y, 



fljv,T(0,O;i,y) oc e- |y|2/2t My) 



dzf N (T-t,z\y) 







j dzsgn 


h N (z) 


j dzsgn 


h N (z) 



det 

l<j,k<N 
p -\z\ 2 /2T 



0<z\<---<z N 

e -y 2 ] /2t-(y j -z k f/2(T-t) _ £ -y] /2t~(y 3 +z k ) 2 /2(T-t) 



x det 

l<j,k<N 



exp 



T 



2t(T - t) 



exp 



T 



2t(T-t) 



Then we set oj' k — tzk/T,l < k < N, and regard oj' = (uj[ , • • • , u>' N ) as the nonnegative eigenvalues of the BdG 
Hamiltonian H' = {H' jk ) S M BdG (2N;R). By (U, dz oc dw' oc d«7|Mw')l' and we have 



?7V,T(0,O;t,y) 
oc My) / dU 



I e -TW\ 2 /2e- det 

h N {oj') l<],k<N 



exp 



r 



2t(T - i) 



exp 



T 



(y 3 + u' k ? 



2t(T - t) 



• (41) 



The result recently reported by Nagao |2J| will give a version of Harish-Chandra (Itzykson-Zuber) integral formula 
in the present case, 



dU exp ( --^tr^W - Ti! 



1 



oc 



det 



h N (u)h N (u') i<i,fe<w 



e -(u]j-uj' k ) 2 /2o 2 _ £ -(wj+Wfc) 2 /2a 



where the integral is taken over the unitary matrices U such that ill £ Sp(2iV; C), and TL and Ti' are hermitian matrices 
in .M BdG (2iV; C) and A4 BdG (2N;TVj having the nonnegative eigenvalues u = (wi, ■ • • , wjv) and uj' — (uj[, ■ ■ ■ ,lu' n ), 
respectively. Application of this identity to (|41|) gives 

3jv,t(0, 0; i, y) oc My) 2 J dU J dW exp f- -l^triH') 2 ^ exp (-S-tr(U^YU - Ti') 



r 



with a 2 = t(l - t/T) and (cr') 2 = i 2 /T, where Y = 
dia,g(y 1 ,---,y N ,-yx,---,-y N ), yj > 0, 1 < j < TV. This 
can be regarded as a BdG version of the two-matrix 
model studied in jlOj for the vicious walker model with- 
out wall. Since it is a convolution of two Gaussian dis- 
tributions, we will arrive at the equality 

g N , T (0, 0; t, y) oc My) 2 / dUJ2 N , T (t, U^YU), (42) 



where, for H G X BdG (2iV; C), 
Mjv,t(*j Ti) 

■ e -{al) 2 /2ol e ^) 2 /2al ^(b^) 2 /2a 2 



N 

n 



y^al v/2^ 2 " v/2^ 2 



n 

l<j<k<N 



,-(«f») a 7»S e -(<4) 2 /- 2 
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n 

l<j<k<N 



(43) 



with 



a' + (a'y = t, <7 ( = a 2 =t[l-- 



Now the transition from the class C to the class CI is 
explicitly represented. The variance <r^ increases linearly 
in t, but of increases in time only up to time t = T/2 and 
then decreases in time. As t — > T, of — > making the 
imaginary parts of matrix elements zeros with probability 
one, and the symmetry class is changed. 



V. MOMENTS OF VICIOUS WALKERS WITH 
A WALL 



Note that 

tr(H )= %h32%hj3 ' ' ' "^■hk-X32k^02k0l1 

where the sum is taken over all N 2k combinations of in- 
dices ji, j 2 , ■■■,32k- 

Since l|43[) is a product of independent Gaussian 
integration-kernels, we can apply the Wick formula with 
the variances 



((4) 2 > = t (1 + Sit) > ((o * )2> = i (1 ~ Sje) ' 

((bfe) 2 ) = %1 + 5 jt ), {{b) t f) = a l{\ + S jt ), 
for 1 < j < £ < N, where 5ji is Kronecker's delta. These 
relations are rewritten as 



In this section, we study the moments of positions of 
vicious walkers with a wall in order to characterize the 
transition of distribution, as reported in [13J for the vi- 
cious walkers without wall. The n-th moment of the po- 
sitions of walkers is defined as 



N 



m N 



T (t,n) = (J2^ 



\3=l I t 
N 



= dx J2x]g N , T (0,0;t,x) 

J0<X 1 < - <X N 

for n = 1,2, • • •, where Xj denotes the position of j-th 
walker. 



A. Wick's formula 

By (|39|) and the equality (|42|) . if n — 2k is even, k 
1, 2, • • •, we have the equality 



where 



m N , T {t,2k) = -(trH 2k ), 



(tyH 2k )= J trH 2k jlN,T(t,H)dH. 



(44) 



(a,j£a mn ) = (aji(a T ) nm } — —(5j n Si m + jSj m di n ), 



(bjib 



(bje(b^) mn ) — —{Sj n Se m +Sj m S£ n ), 



for 1 < j, t,m,n < N, where 



, t(2T - 1) t 
c = --^ — -, 7 



T 



2T -t 



Define 



+ $j+N mSe-N n + ^j-Nm^i+Nn- (45) 

Then we have the variance of the BdG-type Hamiltonian 
in our time-dependent ensemble as 

2 

(TijiHmn) = —^5 N {j,l;m,n) + jS N (j,£;n,m)^ (46) 

for 1 < j,£,m,n < 2N. The Wick formula for 101} is 
thus 



m N,T 



(t,2k) — ^ (^J w (l)Jw(l) + l^ijr(a)i«(2)+l)('^jV(3)i w (3)+l'^2ir(4)iir(4) + l) 

(^iir(2fc-l)J-ir(2fc-l) + l ^jw(2k) J»r(2fc) + 1 ) ! 



(47) 



with the identification j 2 k+i = ji, where the first sum is taken over all N combinations of indices j\,j 2 , ■ ■ ■ ,j 2 k, 
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FIG. 3: Ten kinds of ribbons with and without arrows. 
Weights are also listed. 




FIG. 4: An example of allowed ways of putting arrows on 
ribbons. The O-vertices are marked by circles. 



and the second one over the set of permutations S^k of 
{1, 2, • • • , 2k} with the restriction 

R: tt(1) <tt(3) < ••• < ir(2k- 1), 
7r(2j - 1) < w(2j), 1 < j < k. 

The total number of the terms in the second summation 
is (2k- 1)!!. 



B. Mobius graph expansion 

There are ten terms in the variance 146|) with l|45|l . We 
will represent each of them by a pair of lines as shown 
in Fig. 3, by expressing Kronecker's delta 5j n by a line 
without arrow connecting j and n, <5j+jvm by a line with 
arrow in the direction from j to to, and 6j-N m by a line 
with arrow in the direction from m to j , respectively. The 
weights, which should be multiplied to the factor c 2 /2, 
are also listed. We regard these pairs of lines as the hems 
of ribbons. As shown in Fig. 3, in ten kinds of ribbons, 
half of them are twisted and others are untwisted. Two 
kinds of ribbons do not have any arrows on hems, and 
other eight kinds of ribbons have arrows on hems. We 
will call the ribbons having arrows in the same direction 
current ribbons (c- ribbons for short), the ribbons having 
arrows in the opposite directions exchange ribbons (e- 
ribbons) , and those not having any arrows normal ribbons 
(n-ribbons), respectively. 

Inserting 14rj|) into (|47|1 gives a sum of the K — 
(2k - 1)!! x 10 fe terms in the form, m N , T (t, 2k) = (1/2) • 
(c 2 /2) fc J2?=iLi with U = T, jlt ...,j 2k LtVu ■ ■ -,32k)- As 
explained in |l3j . each term Lg(ji, • • • , j2k) is represented 
by a graph which consists of a 2fc-gon with its edges 
3x32 , 32]?, , ■ ■ ■ , 32k ji connected by k ribbons to make k 
pairs (Wick pairs). Such graphs may be called Mobius 
graphs having ribbons with and without arrows. If we 
assume that there are k n n-ribbons, k c c-ribbons and 
k e e-ribbons, and among them ip n n-ribbons, (p c c- 



thc weight of 



x 7 



-<Pe - 



ribbons and ip e e-ribbons are twisted, 
Hh, ■ ■ -,32k) is 7 V " x (- 7 ) fc -<^ x (_!)* 
(_l)fec 7 fe-fe„- ¥ >+2 ¥ >„ ; whcrc k = k n + k c + fc e and ip = 
+ + </?e (the total number of twisted ribbons). 
For each vertex j s ,l < s < 2k, we take the summa- 
tion of the index over 1 < j s < 2N to calculate Li from 
{Li(ji, • • • , 32k)}- Since each ribbon represents a product 
of two Kronecker's deltas in H45|) . any pairs of indices j s 
and j' s connected by a line (a hem of ribbon) should be 
identified, or identified in modulus ±AT, and the free in- 
dices remaining after this "identification" of indices give 
N dependence to Lg. We will find the following rules for 
the dependence, where it should be noted that each 
vertex is the endpoint of two lines (two hems of two rib- 
bons) with or without arrows, (i) If both of the lines 
connected to a vertex have no arrows, then we will call 
such a vertex a 0-vertex. The summation over a free in- 
dex on a 0-vertex gives 2A^. (ii) If at least one of the two 
lines connected to a vertex has an arrow and they are not 
a pair of lines with inward and outward arrows, then the 
summation over a free index on such a vertex gives N. 
(iii) Otherwise, the sum over free index becomes zero. 
Consider the equivalence classes of the Mobius graphs 
with n-, c-, and e-ribbons. If a class is expressed by 
a representative graph, say T, the number of elements 
of the class T (i.e. the number of graphs topologically 
equivalent with T) is denoted by Let ^(r) be the to- 
tal numbers of free indices of T and Vo(r) be the number 
of free indices on O-vertices. Set Q(k) be the collection 
of all such graphs {T}, then we have the combinatorial 
expression for the moments as 



tT (t,2k) = l -( c -\ £ ir^n^r) 



m N > 



ree(fc) 

x , 1 ^k c (r)^k-k n (r)-<p(r)+2 Vn (r)_ 

In 0| the collection Q of topologically distinct Mobius 
graphs, which consist of 2fc-gons and k normal ribbons, 
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was introduced and the fact was used that each graph 
r 6 Q having only untwisted ribbons (resp. having 
some twisted ribbons) defines an orientable (resp. non- 
orientable) surface Sr by a map [3, E3 to derive the 
1/N 2 (resp.l/A) expansion. The number of distinct 
orientable (resp. non-orientable) surfaces with genus g 
obtained from the graphs without any twisted ribbons 
(resp. with m twisted ribbons) in Q was denoted by 
e g (k) |3g (resp. £ g , m {k) ,131). We notice that the 
Mobius graphs in Q{k) introduced here are obtained by 
putting arrows on hems of some of the normal ribbons in 
r € Q{k). There only allowed ways to put them are such 
that there is no pairs of lines with inward and outward 
arrows connected to a vertex (see Fig. 4). We now in- 
troduce the following multiplicative factors to e g (k) and 
Sg,m(k), A% Kn Kc (k) = the number of allowed ways to 
put arrows on the hems of ribbons in a Mobius graph 
without twisted ribbons, which is mapped to a surface 
with genus g, so that the graph becomes to have K n n- 
ribbons, n c c-ribbons, and vq free indices on O-vertices, 
and A 9 ,^ Kr m (k) = the number of allowed ways to put 
arrows on the hems of ribbons in a Mobius graph with 
m twisted ribbons, which is mapped to a surface with 
genus <?, so that the graph becomes to have n n n-ribbons 
in which m„ are twisted, n c c-ribbons, and Vq free indices 
on O-vertices. Here remark that the total number of rib- 
bons is fixed to be k, the number of c-ribbons should be 
k — K n — k c , and that the total number of free vertices 
V(r) is determined by k and g through the relations with 
the Euler characteristics x = V — k + 1 as x = 2 — 2g 
for m = and x = 2 — g for m > 1, respectively. Then 
we have the following 1/N expansion formula for the mo- 
ments, 



cases, 



2 fc [fc/2] 

m N , T (t,2k) = - i^-j £ *M [jp) 



k+l-2g k 

x E EE( 4 ) Kcyo ^^J^" s " 

t?0 = K n —0 K c 



1 (<?_ 

2 V 2 



,N, 

g—l v / m—1 



k+l-g 



v =0 



E E E i-^AlZ^raSk)^ 



K n — k c m T , — 



(48) 



Here we may prove that A 9 ,^ Kc mn (k) — for k — n n — 
m + 2m n < 0, and (|48(l will give a series with nonnegative 
powers of 7. 

As an example, we consider the fourth moment. For 
k = 2, we found Q3 

£ (2) =2, ei(2) = 1, 

ei,i(2)=4, ei, 2 (2) = 1, £ 2 ,i(2) = 2, £2,2(2) = 2. 
As we can confirm easily that the factors A i j ]Q Kn Kc (2) and 
A 9 1 ' m K K m ^ (2) have non-zero values only in the following 



^3,2,0(2) — 1; ^1,2,0(2) — 1j 
< M (2)=4, Aj >0)0 (2)=2, 

^iii,o,o(2) = 2, ^'2,0, i(2) = 1; ^o'o,i,o(2) = 4, 

^2*2, 0,2(2) = 1, ^-0^1,1,0(2) = 2, -^1^2, 0,l(2) = 1, 

A)'o, i,o(2) = 2, ^.o'l, o,i(2) = 2, A)'o,o,o(2) = 2, 
^1,0,1(2) =4, ^,0,2(2) = 1- 

Then (O gives 

c 4 r 

m NjT {t, 4) = j{8A 3 + (6 + 87 + 2 7 2 )A 2 

+ (l + 2 7 + 5 7 2 )iv}. (49) 

By definition we will see A® Kn Kc (k) = 
S Vo k+i$K n kS Ka o, and thus the leading term in lf4*%|l 
for large N is 

m N , T (t, 2k) = ~ (j\ (2N) k+1 e (k) + 0(N k ). 

Since so(k) is the Catalan number, Wigner's semicircle 
law will hold in N — > 00 also in the BdG random matri- 



ces. 



C. Calculation by density function 

By (tT2l with (JJJ and ijTTl) . it is easy to see that 
9n,t(0, 0; t, x) is symmetric in x±, ■ ■ ■ , xn- Then we can 
define the density function as 



P(t, x) 



(N-l) 



— — / Y\dxjg N ,T(0,0;t,x) 



and the even moments l|44f) are calculated by it as 



m N , T {t,2k)= I x 2h p(t,x)dx. 
Jo 



(50) 



Let L^ a \z) be the Laguerre polynomials with a param- 



eter a defined as 



L[ a) (z) 



e z 



{e- z z ]+a ) 



Using them with a = 1/2, we can define the monic poly- 
nomials Cj(z) = (—iy j\L^ 2 \z), which satisfy the or- 
thogonality Jq 00 z 1 / 2 e~ z Cj(z)Ct(z)dz = hjSj£ with hj = 
T(j + 3/2)j!. Quite recently Nagao gave the general ex- 
pressions of dynamical correlations for the present system 
using Cj(z) |24| . From his result, we can read the density 
function p(i, x) as 
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3=0 



oo JV-1 I 



where 



/3j'2<? 



(-i) j+1 r(i- 2^-1/2) j\ 



27 



i\ r(2j-^-i/2) 



4/ (2j-£+l)W. ' 



, _ (-1)' jj [(J 2 ^ 1)/2] r(j" - 21 2m - 3/2) 

2^7? (j-2£)\ (2£)V h2£+1 2^(2£+l)\ ^ (j - 21 - 2m - 1)! ' 



Let Hj(z) be the j-th Hermite polynomial, Hj(z) = (— l)- 7 e z (d/dz) J e z . If we use the following equalities 37, 38], 
LC- 1 ^) = Lf\z) - L^z), zL ( ; +1 \z) = (j+a + l)L^(z) - (j + l)L%{z), 



E 



-^V{m + a+l) 

■m— v 7 



r(j + a + 1) 

' ' 2)(Z) = W^'^' ^4 1/2) ( 2 ) = ^T7T^i + i(Vi) forz>0, 



we will have the following expression for the density function, 



p(t, x) 



N 



1 



247V-3(7v_ l)!r(JV+ 1/2) c 



-(*/c) 2 



irrr f I \ \ 2 (2JV-l)(JV-l)c 



7V 



N —1 c 

H 2 n (x/ c)H 2 n-z{x/c) — H 2 n (x/c)H 2 n-i(x/ c) 

x 2N x 



1 oo N—X i 



j=N £=0 m=0 



2 2m {2j + 1)!0F 



>/ c ) 2 tf 2i+1 (z/ C )if 2m+1 (x/c). 



(51) 



Substituting (|51|l into (|50|l and using the integral formulae of Hermite polynomials used in 0] and the relation 
Y^i=m Pniotim = S nm for n > to |24| , we will arrive at the expression 



m N , T (t,2k) = (^j N(2k- l)!!g (T-/) 2 ^ 

— J EE E 7m+1 E Pm-i+NnOtnN-l-1 
' j=Q 1=0 m=i n=N 



k - 1\ ( k-j-1 



N - 1 



j ) \2(k-j + l) 2N-k+j 



{2k)\{2N -2£- l)\2 m - j+1 



' + 2N - 2£ + m - k)\(k - j-m- l)\(k - j + m + 1)! 

I 



(52) 



As a matter of course, when we set fc = 2, l|52|l gives 
the fourth moment l|49|l . The large- N behavior discussed 
below jUDl is obtained also from (JS^J. 

VI. CONCLUDING REMARKS 

In the present paper we have considered the solv- 
ability of one-dimensional vicious walker models, which 
are generalizations of the model studied in earlier pa- 
pers 0, 0| . We are interested in these systems as 



the noncquilibrium statistical models, since they provide 
in general spatially and temporally inhomogeneous sys- 
tems. The temporally inhomogeneous one-particle sys- 
tems have been extensively studied by Yor by construct- 
ing them from Brownian motions and Bessel processes 
[l7J . Our present work may be regarded as a new attempt 
to construct many-particle systems in one-dimension us- 
ing such temporally inhomogeneous processes as elemen- 
tary processes, so that we can discuss also the spatial 
inhomogeneity. 

We have imposed the noncolliding condition between 
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particles, since we have thought them as vicious walk- 
ers. This condition introduces the effective long-ranged 
repulsive interactions among particles and it may make 
the models be nontrivial many-particle systems. Our 
strategy to analyze the models is to map these inter- 
acting particle systems to multimatrix models defined 
in the appropriate matrix spaces, which have in general 
higher dimensions than the original phase spaces of par- 
ticle systems. The particle positions are expressed by the 
statistics of ei genv alues [2Sj or the "radial coordinates" 
CI 111 HJ of random matrices. 

We have proved the following equivalences in probabil- 
ity distributions. 

(i) The noncolliding system of Brownian motions with 
a wall, defined as the continuum limit of the vicious 
walker model with a wall, and the noncolliding sys- 
tem of Brownian meanders, which are constructed 
from the three-dimensional Bessel process. 

(ii) The noncolliding N generalized meanders con- 
structed from the d-dimensional Bessel processes 
with d = 2, 4, 6, • • • and the "radial coordinate" of 
the {N+ (d— 2)/2}xJV rectangular random matri- 
ces in the Gaussian ensemble with time-dependent 
variances (a two-matrix model of the chiral GUE 
and chiral GOE.) 

(iii) The noncolliding N Brownian meanders and the 
nonncgative eigenvalues of 27V x 2N BdG-type ran- 
dom Hamiltonians in the Gaussian ensemble with 
time-dependent variances (a two-matrix model of 
the BdG-type random matrices in the classes C and 
CI of Altland and Zirnbauer). 

We have discussed the noncolliding system of general- 
ized meanders made from the general odd-dimensional 
Bessel processes. So a natural question is what is the 
corresponding random matrix theory for the noncollid- 
ing system of generalized meanders associated with the 
d-dimensional Bessel processes with d = 5, 7, 9, • • •. As 
far as we know, it is an open problem. 

We have claimed that the above equivalence between 
the interacting particle systems and the Gaussian mul- 
timatrix models implies the solvability of the systems. 
This statement may be true, but to obtain exact expres- 
sions of general correlation functions for the multimatrix 
models |39l E3 is far from trivial and one have to use a 



series of techniques developed in the random matrix the- 
ory. Exact expressions of general dynamical correlation 
functions enable us to discuss the infinite particle limit 
N —> oo of the nonequilibrium system as reported in 
[l2l l4lj for the vicious walker model (non colliding Brow- 
nian motions) without wall, and in |24j for the system 
with a wall. We expect that the infinite particle limits 
of dynamical correlations for the noncolliding general- 
ized meanders can be generally evaluated by following 
the strategy employed in 0> 03 • 

As a simplest case of the correlation functions, the den- 
sity function can be determined. As reported in |l3[ . cal- 
culation of the moments of the positions of vicious walk- 
ers is related with an enumeration problem of orientable 
and non-orientable surfaces with a fixed number of genus, 
which are obtained from the Mobi us grap hs with a fixed 
number of twisted ribbons by map |34ll35j| . In the present 
paper, we showed that a new graphical problem arises 
from the vicious walker model with wall restriction; an 
enumeration problem of the ways of assigning arrows on 
the ribbons of Mobius graphs following some rules. It 
should be noted that, roughly speaking, there were two 
kinds of ribbons in our expansion formula, with and with- 
out arrows. These ribbons are thinned into lines and they 
are drawn on surfaces. Then this enumeration problem 
may provide, if we consider the lines as "world lines" of 
particles, the statistical mechanics of composite particles 
on random surfaces. 

Recently a variety of problems associated with the con- 
ditional random walks /Brownian motions have been pro- 
posed and intensively studied in statistical physics, e.g . 
first passage problem , "lion-lamb" problem [43 , 
diffusion particle systems with mobile traps ^(| , families 
of vicious walkers Q, leader and laggard problem |47) . 
system of stochastic Loewner evolutions |4jj, friendly 
walker models [M |H HJ H El and so on. Solvabil- 
ity and unsolvability of these models will be important 
topics in statistical physics far from equilibrium. 
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